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POTENTIAL 
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Abstract. We give two algorithms that allow to get arbitrary precision asymp- 
totics for the harmonic potential of a random walk. 



1. Introduction 

The discrete harmonic potential a of a random walk R, starting from 0, on a 
lattice Z can be most easily defined using 

oo 

a(z):=^P„(0)-P„(z) (1) 

Tl = 

where Pn{z) is the probability of R to be at z at the nth step. It is easy to conclude 
from Q that Aa = (5{o} where A is the discrete Laplacian for the walk R. Less 
obvious are the facts that a is positive, and that these two properties determine 
a up to the addition of a positive constant. "Positive" may be replaced by "with 
sub-linear growth" in dimension greater than 1 . 

The discrete harmonic potential is an interesting quantity, strongly related to 
the discrete Green function g{x, y): on the one hand, a{x) = g{x, x) where g is the 
Green function corresponding to the set {0} {576} P11.6, page 118]. On the other 
hand, g for arbitrary sets can be calculated if a is known — when the set is finite 
there is an explicit formula, IS76I T14.2, page 143]. 

Thus it becomes interesting to calculate or estimate a. Except in dimension 1 
where a{x) = \x\, the problem is far from trivial. Take as an example the harmonic 
potential of the regular random walk on which has a logarithmic nature. The 
estimate of a most common in the literature is 

a(z) = -log|z| + A + 0(|zr2) (2) 

where A is some number, which can be expressed using Euler's 7 constant, A = 
■2-7 + i log 8. The earliest proof of ^ we are aware of is in Stohr LS49J . A more 
accessible proof of a weaker result, giving an error estimate of o(l), can be foimd 
in Spitzer IS761 P12.3, page 124]. Spitzer's proof (which is not specific to the lattice 
Z^), like Stohr 's, relies on the fact that the Fourier transform P„ has an explicit 
formula, which can be summed to give a "pseudo-Fourier" representation of a as 
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the expHcit integral 



ai.)^J-^d, ; (3) 

indeed, this is how he proves that the sum in (0 converges in the first place. Fukai 
and Uchiyama used this technique to find the next order approximation — in the 

case of the simple random walk on it is — g^^fp — and to show that an asymp- 
totic polynomial expansion exists. See IFU96I for the two dimensional case and 
lU98l for the case d>3. 

The purpose of this note is to give two alternative approaches to the computa- 
tion of a, which allow to get asymptotics of arbitrary precision. For example, here 
are the first few terms for the regular random walk on Z^ : 

2 / 
a{z) = — log l^l + A — Rc 



TT °' ' ' \67r|zp 



+ 



207r|z|8 247r|zP 
51z^ 352^2 



567r|z|i4 367r|z|i8 
217z*^ 45zi2 1925zi6 



1607r|z|i6 47r|z|20 1927r|2|24 

-10\ 



+o(izr"')j • (4) 

The first approach is a brute-force one which starts from O: P„ may be repre- 
sented as a sum of multinom coefficients, and plugging in Stirling's expansion, 
and some elementary algebra gives the representation JIJ. Unfortunately, the al- 
gebra involved is too long. In effect, for all but the first coefficient (the —^^^), 
it is simply not practical to do the calculations by hand. We demonstrate this ap- 
proach in two ways: we use it to prove that a pol3momial approximation exists, 
reproducing the results of IF U 961 IU98I in the case the random walk is bounded 
(Fukai and Uchiyama show these results under the weaker assumption that each 
step of the random walk has K + d + e moments, where K is the required precision 
in the expansion, d is the dimension and e > 0). See theorem^on page|5j and we 
use it to get some explicit constants and high order expansions for other walks. 
See sectional 

The second approach is to approximate the discrete Laplacian A using an ap- 
propriate differential operator. For example, for the regular random walk on Z^ 
we have, from the Taylor expansion, 

(A/)(z) - i(/(z + l) + /(z-l) + /(z + z) + /(z-z))-/(z) 

(5) 







1 fd\f 






dy'^ J 


^ 48 \dx^ ^ 


dy^J 



In fact, A = i(cosh^ + cosh^) — I, but we will not use this representation. 
Given that the expansion of a has an appropriate polynomial form, the expansion 
l|5j allows to get many equations on the coefficients in (4j. These equations are 
of course insufficient — the equation Aa = S^o} does not determine a uniquely, 
and the condition that a is positive or slowly increasing is hard to encode into the 
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differential equations. However these differential equations allow to prove many 
of the apparent properties of most notably that the coefficient of z^\z\^^ for 
/ < ^fc is always zero. This is theorem|2 Unfortunately, it is still not clear why all 
the coefficients are negative. 

Returning to the case of the regular random walk on I?, it can be seen that the 
equation Aa = combined with the -^-rotational symmetry gives that the val- 
ues of a are uniquely defined by their values on one diagonal, say {m + i'tiIJ^^q, 
and vice versa, any arbitrary sequence on the diagonal can be extended to a sym- 
metric function / with A/ = J{o} using a very simple recursion. On the other 
hand, it turns out that the values of a on the diagonal can be calculated explicitly 
using the integral (3| to get 

a(.n + ^m) = l(l + i + l + ... + ^) (6) 

IS76I chapter 15]. Indeed we have just sketched an algorithm for calculating any 
value of a(z): this is the McCrea- Whipple algorithm IMW40 I^. As is now evident, 
all the values are of the form n + -q, n G Z and g G Q. Interestingly, when not 
on the diagonal, n and q increase exponentially, even though the result is only 
logarithmic in size^. As for the coefficients of l|3J, it turns out — perhaps unsur- 
prisingly — that l6j can be used to complete the few coefficients that cannot be 
deduced from Jsj. This is theorem|3l 

There is a third approach to the problem which we shall not discuss at length as 
it seems inferior to both former ones. Roughly it goes as follows (for the case of the 
regular random walk on I?). "Guess" that the solution is approximately ^ log \z\. 
We discretize to 1? in the natural way: define 

/(^) 

Actually, it is not necessary to guess the exact value |, other constants not too 
different also work. A calculation can show that 

^ |(A/-<5{o})WI<l . (7) 
This allows us to write a series of corrections of / as follows: /i = / and then 

fn = fn-1 - (A/„_i - (5{o}) * / 

where * is the convolution on Z^. Since A{g * h) ^ g * Ah we get 

A/„ - 5(0} = (A/„_i - 5{o}) * (A/ - <5{o}) 

and hence l|7} gives that the /„'s converge^ in the norm to an a which satisfies 
Aa = (5{o} and \a\ < C log |z| + C so it is the harmonic potential up to an additive 
constant. Furthermore, the fact that A/ — (5{o} = can be used to derive 




^ Can be found in Spitzer as well: IS76I chapter 15, page 148] 

^ No, this is not a very efficient way to calculate tt, it requires at least operations to get a precision of 
k digits... 

^ The convergence rate is obviously exponential, but one might also define fn = fn-i ~ {A/,i_i — 
(5{Q-j ) * /n-i and get double exponential convergence. 
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similar estimates for a. However, the best approximation we were able to get from 
that method is 

a{z) = I log |z| + A - Rc + 0{\z\~^ log \z\) (8) 

and getting this approximation was not significantly easier than the brute force 
approach. Moreover, this approach did not allow to get actual values for any of 
the constants — only that some constants exist, and their values were derived from 

We have described the contents of this paper except for the last section. That 
section contains the results of some computer-aided simulations and a calculation 
of an explicit constant in the 0{-) in O. 

We wish to thank Greg Lawler for some useful remarks. 

1.1. Standard definitions. A lattice is a discrete additive subgroup of M". The 
dimension of a lattice Z (denoted by dim Z) is the dimension of the linear span of 
Z as a linear subspace of R" . A basis for a d-dimensional lattice is a set {ei , . . . , e^;} 
such that Z = e^L + • • • + e^jZ. Z' is called a sublattice of Z if it is a subgroup of Z, 
and its index is the size oi Z/Z' . The volume of (the cell of) a lattice Z, denoted by 
vol Z, is the volume of M'^/Z (the discreteness of Z allows us to select a measurable 
set of representatives). Alternatively it can be defined as 

r^oo #{B{r,0) n Z) 

where B{r, 0) is a ball of radius r around and \B{r, 0) | is its volume. 

A random walk i? on a lattice Z is a probabilistic process {Rq, Ri, ■ ■ ■}, Ri G Z, 
such that i?o = and Ri — Ri-i has a distribution independent of i, and such 
that for every z Z there exists some n such that P(z = i?„) > 0.* It is bounded if 
Ri~Ri^i is bounded. We will typically confuse the process R with the distribution 
of any step (e.g. with the distribution of i?i) so that we can use notations such as 
Ei? comfortably. The dimension of R is the dimension of the lattice Z, and is 
denoted by dim R. 

The drift of a random walk R is E_R. The random walk is balanced if its drift is 
zero. 

The discrete Laplacian A of a random walk R is an operator on functions on the 
lattice Z defined by 

(A/)(z) = E/(z + i?i)-/(z) . 
Fimctions with A/ = are called (discretely) harmonic. 

The continuous Laplacian on Mr, + • • • + will be denoted by Ac. Ftmc- 

tions with Ac/ = will be called continuously harmonic. 

I |u| I will always refer to the L? norm. Similarly, for a matrix A, \ \A\ \ refers to its 
norm as an operator from to IJ^. 

2. The direct approach 
The aim of this section is to prove the following theorem. 

* This definition is a bit restrictive. For example, a 1-dimensional random walk which always goes to 
the right is not a "random walk on a lattice". 
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Theorem 1. Let a be the harmonic potential of a d-dimensional balanced bounded random 
walk on a lattice Z. Let K > be some integer. Then there exists a constant t, a linear 
map A, a polynomial Q ~ Qk in d variables and an integer L = Lk such that 



a{z) = T^,a{\\Az\\) + + 0{\\z\\-^) (9) 

where 

fJ'diz) := 



\og\\z\\ d = 2 



T and A have expHcit formulas — see the comments after the end of the proof 
(paeelTOlbelow). The expression (5(z)/||A2:||^isA + 0(||Az||^^'') (this is classical, 
of course, but will also follow from the proof below). The number A is of course 
also walk-specific. 

In the case d = 1 the equation Aa = (5{o} deteriorates into a simple recursion 
formula and the theorem is nothing but an exercise (and Q = A). Thus we concen- 
trate on the case d > 1. 

Lemma 1. Let pi, . . . ,pk G ]0, 1[ such that = 1. Let K £ N. Then there exists a 
Q = Qp,K such that 

Bp{n,pin + wi,...,pkn + Wk) = (10) 

(27m)('^ ^i^^^px ■ ■ -pk 
for every wi,. . . ,Wk such that ^ wi — 0. B here is the standard multinom 

n ' 

Bp{n,mi, . . .,mk) = j rPi Pj. , mi -\ \- ruk = n 

nil'. ■ ■ ■ ■ mk'. 

and Q satisfies that n^Q is a polynomial for some L = and Q{n,w) = 1 + o(l) 
uniformly in \ \w\\ < C^Jn log n. 

The proof of theorem ^ is replete with these rational functions Q obeying the 
condition above. Notice that this simply means that 



Q{n,wi, . . .,wk) = ^ 



c.i7n^Wi^ ■ ■ ■ wl. 



such that = 1 and otherwise c-j =/= only if -L < z < and < ji + • • • + jk < 
— 2i. It turns out (and this can be deduced from the proofs below with some care) 
that Lk = 4A' is always sufficient. However we will have no use for this fact. 

Proof of lemma^ Consider Stirling's series, 

n\ = — V2nn-exp {Qi{n) + 0{n-^^)) 
where is some polynomial with Qi{n) = o(l) as n ^ oo. This gives 

1 A 1 



Bpin,pini +wi,... ,pknk + Wk) = ——rj-^y^-== J] 



{2nn)ik~i)/2^p^...p,lL r 



Pin 



exp [ - ^{pxTi + Wt)log ( 1 + — + Qi{n) - ^ Qi{pin + w,) + 0{n ^) 



i=l 
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Assume for the moment that | | < C\/n log n where C will be fixed later. In this 
case we can insert monomials such as \k\\ into 0(iir ). Thus, each of the factors 
, ^ equals Qiin, wi) + 0{n^'^) with n'^^Q2 some polynomial with Q2 = I + 

V 

0(1) {Q2 depends on pi, of course). The term inside the exponent evaluates to 

2 

-Y.-^^ + Q3{n,w) + 0{n^^) with n'^'^Qs a polynomial satisfying Q3 = o(l) 
uniformly in lluijl < Cy/n log n as n ^ 00. This proves the case ||in|| < C^/nTogn. 
In the case ||w|| > C\/?i log n both sides of ITUt are 0{n^^) which follows easily 
from the fact that in this case, for C sufficiently large 

,,2 



exp <cxp(-Clogn)<Cn"'f^ . □ 

Lemma 2. Let v GM.'^andAe GL(R'^). Then 

^exp(-l||^(z + .)|p) =M!!^! + 0(e— ) 

where the O(-) is uniform in v. (\A\ here is the determinant of A). 
Furthermore, if P is any polynomial and y <eW^ then 

^ cxp('-i||A(z + «)|p')p(z)=n'^/2Q(„) + 0(e— |H|dcgP) (11) 

where Q is a polynomial. The coefficients ofQ are fixed polynomial functions (depending 
on A) of the coefficients of P and of v. Q is independent ofy and the O(-) is uniform in v 
and y. 



The use of the parameter n in the formulation of the lemma is a bit artificial. 
Later, when we shall use the lemma, n will be an integer number (the number of 
steps), but this fact is not needed for the formulation or proof of the lemma. 

Proof. We use the c?-dimensional Poisson summation formula, J2 ./(") /(27i'«) 

(the Fourier transform being defined by /(^) ~ J /(a;)e*^^''>Ma;) for /(z) = cxp (— i| + w)| {z+ 

vY where r=: (ri, . . . , r^) G {0, 1, . . . j'' and {z + vY is short for (zi +viY^ ■ ■ ■ {zd + 

VdY''- A simple calculation shows that 



/(^)=,-.(..),-.|_(![!^exp(-in {A-r^ 



(12) 



where r = ri + - ■ ■ + rd, and as usual = -^rn- ■ ■ ■ -rrnr- This means that as n ^ 00 
we have 

(note that \A\~^ < WA'^W^). The exponent inside the O(-) kills of course all the 
other factors. The case r = immediately gives the first part of the lemma. For the 
second part, put ^ = in 1 1121 and get that /(O) 7^ only when r is even, and that 
/(O) = n'^/'^Qp{n) where Qp depends only on A. Writing P = J2r'^r{z + I'Y we 
get Q ~ J2r'^rQr which concludes the lemma: the fact that cp are polynomial in 
the coefficients of P and in v gives the same for the coefficients of Q; the fact that 
cp = OdlvW^"^^) gives the same in the O(-) in The result does not depend 
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on y (except in the O(-) error) since f{n) does not depend on w — in effect the 
dependence on v appears only via the dependence oi cponv. □ 

Lemma 3. Let Rbea d-dimensional balanced bounded random walk on the lattice Z. Let 
K G N. Then there exists a sublattice Y c Z x Z, a constant t', a quadratic function B 
and a Q such that 

jT'n-^/\-B(^y^^Qin,z) + 0{n-'<) {n,z)eY 
10 (n,z) 

with n^Qa polynomial in nand in the coordinates of z, and Q{n,z) = l + o(l) uniformly 
in \ \z\ \ < C^Jn log n. 

Assume for simplicity that Z c M''. A simple comparison of the expectations 
of the left and right side of l IT^ shows that B{z) ~ ^{M^^z, z), where M is the 
correlation matrix of R, Mij = K{R, ei){R, ef), where e; are unit vectors. Notice 
that the expression (Af^^z, z) is independent of the coordinate system chosen. 

The lattice Y is clearly d+l dimensional since otherwise we would have dim R < 
d. Thus Y O Z is rf-dimensional and vol Y D Z is finite. This, with the requirement 
^ P„ = 1 allows to calculate r' and get 

r' = (14) 

We will not be using these equalities in the proof, though. 

Proof. By applying a linear transformation we may assume i? is a random walk on 
Z^. Assume R moves from z to z + Vi with probability pi for i ^ 1, . . . , N + 1. Let 



W 



|(n, wi, . . . , wn+i) ■ Pin + Wi G Z Vz and Wi = o| 



which is an {N + l)-dimensional lattice. For an n e N let hi be the number of times 
(between and n — 1) that R moved from z to z + Vi. We use lemma^with some 
K2 which will be fixed later. This gives us 



"{ki = npi + w,) = Cn cxp ^-^B2{w)J Q2{n, Wi,. .. , Wn) + 0(n 



K2 



) (15) 



where B2 is a strictly positive quadratic form, C is some constant and n Q2 is a 
polynomial with Q2 = 1 + o(l) imiformly for \ \w\\ < Cy/n log n. Now, for a point 
z eW^ we can write 

P(i?„ = z)= Hh = np, + w,) (16) 

W{n,z) := |?« : (n^w) EW A "^^WiVi — z| 

Since W{n, z) are intersections of a lattice with a translated subspace, we get that 
W{n, z) = for (n, z) outside some lattice Y , and for (n, z) eY we get that W{n, z) 
are translations of one fixed lattice X. Since we know Y is d + 1 dimensional (see 
the comment just after the statement of this lemma) we get that X is {N — d)- 
dimensional. The calculation of the sum in will be done of course with lemma 
|2but we need some preparations first. By applying a linear transformation to the 
A^-dimensional space of w's we may assume that 

W{n, z) = ziei 4 h z<je<j + e<j+iZ H h e^Z + y„ 
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where ?/„ e spa,n{ed+i, ■ ■ ■ , bn}- Combining lITSl , IT6l and the linear transforma- 
tion we get 

P(i?„ = 2)= Cn'^/^exp(-^B3iw)]Q3{n,w)+0{n-''') (17) 

where K3 = — + rf. This last error estimate follows easily from the following 
obvious estimates: 

From now on it would be easier to replace i?3(ui) with for some A e 

GL(R^) which is possible since is also strictly positive. We denote x = (0, . . . , 0, 
Wd+i, . . . , wn) so that w ~ X + z and we can write 

\\Aw\\^ = \\Ax\\^ + \\Az\\^ + 2(x,A*Az) . 

Since A* A is strictly positive we get that its lower right N — d minor is also strictly 
positive and thus invertible. Therefore there exists some v — (0, . . . , 0, Vd+i, • ■ ■ , f jv) 
such that 

{A*Av)i = {A*Az)^, i = d+l,...,N . 
This of course implies {x, A* Av) = {x, A* Az) and then 

\\Aw\f ^\\Az\f ~\\Av\\^ + \\A{:x^v)\\^ . (18) 

The importance of ilSt is that its left part, IIAzjp — is constant for w e 

W(n^ 2), while exp (— i||A(a; + w)| P) can be summed by using lemma|2on the last 
N — d coordinates. Notice that it depends on the z's only via v, which is linear in 
z. In other words, we can now define our quadratic form B: 

B{z) -.^WAzW-WAviz)]]' 

and apply lemma 12 We get 

|T7I) Cn-^/2g-B(.)/„ cxp(^-\\A{x + v)\\AQ3{n,z,x)+0{n-'^') 

= Cn-'^/2e-^(^)/"Q(n, z) + 0{n-'^' ) + 0{e-^'' \ \v\ I'^^sQa) 

We may already remark that picking A'2 ~ K + N — d will give us a precision 
of 0{n^^) as required. Obviously 0(e^'^"||i'||'^°s'5^) = 0{n^^) for any K as v is 
linear in z and | |z| | < C^/n]ogn. To say something about Q, we need to consider 
Qs as a polynomial in x with coefficients which are polynomials in z and negative 
powers in n, i.e. 

Q3 = y^^cp{n, z)x'' 

r 

where n^c^ are polynomials. Lemma |2 tells us that n^Q is a polynomial in n 
and that its coefficients are polynomial in the coefficients of Q3 and in v (which is 
linear in z). Thus we get that n^Q is polynomial in both n and z. This concludes 
the lemma: the only claim we haven't shown is that Q = 1 + o(l) which can be 
deduced, for example, from the central limit theorem^. □ 



^ Actually, it is also possible to deduce it from the arguments above with a little care as to the exact 
dependence between P, Q and v in lemma|5| 
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Remark. For specific walks it is possible to prove lemma [S] using much simpler 
techniques. For example, for the simple random walk on 7?, we have 

Pn{x + iy) = A-'^^ ' ^' ^' 



fe=0 



{^\kj\{k~x)l2) \{n-k-y)/2 



where we make the notational convention that (^) = if a is not an integer. This 
formula allows to get lemma|3]from lemma^more-or-less directly, with no need to 
go through lemma|2] This works for other walks as well. However it does not seem 
to work for the very interesting case of the random 6-walk on the triangular lattice, 
a walk which can be considered no less natural than the random walk on 1? . We 
mean here a two dimensional random walk that goes at each step with probability 
i from z to one of the points z + a;% i = 0, . . . , 5 where uj = ^ + t-^ = ^/l. 

Lemma 4. Let </ e N and 1 < j < q- For every s > 1 and for r oo zve have the 
superpolynomial estimate 

n=j (?) 
n>j 

while for s ~ live have 



"=j (9) 

n>j 



n 



Proof. Again we use the (one dimensional) Poisson summation formula, this time 
for the function 

f a;-"e-i/^ a; > 

^("-■^ = 1 n n ■ 

10 x < 

Simple Fourier operations show that 

Since / e we get that /(O = 0{C^) for all K and hence F = ^/(O) + 

0(7'^^^). This concludes the case s > 1. For the case .s = 1 we notice that the O(-) 
in \19\ is uniform in s S (1, 2] since it depends only on the norm of the /v 'th 
derivative of /. Thus we need only to calculate 

where C is Riemann's function. Writing F = i + C + o(l) and Q{s) = + C + o(l) 
and a little use of I'Hopital rule will give log r + C and conclude the lemma. □ 

Remark. The exact asymptotics of the error are harder to figure out. This is of 
course a question about the asymptotics of /. A much better estimate than /(^) < 

can be had by changing the path of integration of x^'^e^'^/^^^^^ from [0, 00) 
to [0, 1 + z] U [1 + I, 00 + i) or to its conjugate (depending on the sign of ^). On 
the segment [0, ^"^^^(1 + i)] we have e^^/^' = 0{e^^^). On the reminder we have 

= 0(e-'=^). So we get f{C) = 0(e-=^^). 
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In the other direction, there are known precise results connecting the "depth" 
of the zero of / at 0, i.e. the rate of convergence of / to zero, with the property that 
/ is in some weighted W space. See theorem 1.1 in IMS02I . We omit the details. 
This shows, for example, that /(^) ^ 0(e~^^'^^') for any e > 0. 

Proof of theorem^ Let B, r' and Y be the quadratic form, constant and lattice given 
for R and A' + 1 by lemma [S] Denoting 

~a{z):= r'n-"'^ ^Pr^iz) 

we see that a[z) — a{z) is a constant (indeed, the precaution of adding r'nT'^/'^ 
is only necessary in the case = 2 in order to make the sum converge). The 
fact that R is boimded shows that P(i?„ = z) is zero for ?i < c||z|| and of course 
E„<c||z||e"^^^^/" = 0(e-=ll^ll). Therefore we can use O to get 

~a[z)= r'n-^/'(e-''(^)/"g2(n,z)-l)+ ^ O(n-^-i) + ©(e-^H^H) 

n:(n,z)eY n>c||2|| 

Lemma|3]allows to calculate the left sum, while the right one is obviously O (| | z 1 1 ~ ^ ) 
This concludes the proof of theorem^ □ 

We see now that the quadratic form of theoremQ]is the same as that of lemma|3l 
i.e. B{z) ~ ^{M^^z, z) where M is the correlation matrix of R. Thus A ~ M^^l"^ 
— of course, this is not uniquely defined, but the quadratic form 1 1 Az| p is. As for 
T, the factor volF n Z of lemma |3] with the term i of lemma |4] gives volZ (this is 
easy to see) and we get 

volZ f2 
(27r)'^/2|detAf|i/2 ' jr(id-l) d^l 

(the factor 2 in dimension 2 comes from the fact that we formulated /X2 = log | |z| | 
rather than log | |z| P). 



3. The differential approach 

We shall demonstrate the differential approach by proving a few results that 
would seem quite difficult using the direct approach only. 

For simplicity of notation we shall always assume that the correlation matrix of 
R is constant. Thus we call a random walk spherical if it is boimded, balanced, and 

E{R, Ci) {R, Cj) = C6{i, j) (again, e; are imit vectors). 

Theorem 2. Let Rhea 2-dimensional spherical random walk on C. Then 

a{z) ^ alog\z\ + X + Re(j2<^ik^) +Oi\z\-'') (20) 

and aik ~ whenever I < ^k. 

The sum here is over a finite collection of / > and fc > with the condition 

-1 > / - fc > -K. 
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Proof. The complex representation <20t is nothing but a restatement of theorem 
writing .t = ^ (z + z) and y = ^{z — z) we get a representation in the form 



E 



Z Z' 



but of course if both i 7^ and j 7^ then we can cancel one of them with the | z | 
factor. The fact that the end result is real means that aojfc = ajofc and we get <20t . 



qk 



Thus what we need to prove is the fact that aik ~ whenever ^ < |- We change 
the representation again, writing 

aiz) = plogizz) + Pkiz^'z' + 0{\z\-'^) (21) 

k+l>-K 

and now we have to show that (3ki = whenever both k and / are negative (of 
course, (3ki = when k + I > 0). Now, the first thing to notice is that we may 
differentiate illl formally. If R moves to Vi with probability pi then 

(A/)(2) = j;pi/(2 + «,)-/(«) = 

n— 1 i \ i 

where we used the one-dimensional Taylor expansion for the functions f{z + tvi) 
for each i (the constant in the 0( ) above is of course simply 1). Denote the nth 
differential operator in the sum by £>„ and denote the maximum inside the 0( ) 
by A/if . Clearly, 

Dnz'^z^ =Oi\z\''+^-") M„z'=z' = 0(|z|'^-+'-") . (22) 
Now, Aa = with (21) and J23 show that 

( (/31og^^+ E Pkiz'A ^Oi\z\-'') . (23) 

\n=l ) \ k+l>-K J 

This is exactly what we mean when we say that J21t can be differentiated formally. 

Now to the actual calculation. Notice that Di = (because R is balanced) and 
that D2 is, up to multiplication with a constant, the continuous Laplacian, because 
R is spherical. Therefore 

D2{z^z^) = 2-iklz''-^z^-^ (24) 

for some constant 7 7^ 0. Assume that we know I3ki for all fc + Z > — m, and that 
[3ki = if both are negative. We write 

\n = l J \ k+l>-m / 

we get that Sm = if both k and I are negative. Moreover, J23> for K > m + 2 
gives equations for f3ki with k + I = — m, namely (3ki = 2^^k~i.i~i- This shows 
inductively that l3ki = when they are both negative, and the theorem is proven. 

□ 

The case / = in theorem |2 has some potential uses for precise estimates of 
hitting probabilities, so it seems worth of special mention: 
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Corollary. For every spherical random walk on R^, the high-order expansion of the har- 
monic potential a, considered as a function of a continuous variable, satisfies the super- 
polynomial estimate 

/ a{Re''^)dt^A\ogR + B + 0{R-"^) VA'. 
Jq 

Remark. All the above discussion has analogues in dimension > 2, using the no- 
tions of spherical harmonics. Spherical harmonics are functions on S"^^^ which are 
eigenvalues of the angular part of the continuous Laplacian (which is the angular 
momentum squared operator). It turns out that each is a restriction to 5*^"^ of a 
polynomial on R'* and we shall use this representation. For a general introduction 
to the topic, see ISW711 sect. 4.2]. 

• In two dimensions, a factor -P/|z|" (where P is a homogeneous polyno- 
mial) may be decomposed into a sum of the terms z'^'z', k + I ^ dcg P ~ a. 
The high dimension analogue is a unique decomposition into a sum of 
terms Q/Hzjl'' where Q is a homogeneous (continuously) harmonic poly- 
nomial, deg Q — P = dcg P ~ a. 

• The analogue of the fact that only and are harmonic is the fact that a 
term (5/||z||'' is harmonic if and only if /3 = 2degQ + d — 2. This follows 
from the following general formula: 

. ( Q \ ||^|pAcQ+(/3^-/3(2degQ + d-2))Q 

and from the harmonicity of Q (remember that Ac is the continuous Lapla- 
cian). \25\ also gives the analogue of ll24> . 

• The analogue of "k and I are both negative" is " (3 < 2 deg Q + d — 2", or. 



n 



equivalently, "Ag>(Q/| |z| j*^) is never zero, for any 
The analogue of theorem|2]is the following: if R is a spherical random walk 
and 



then 



We omit the proof, which is a similar induction. 

The corollary to theorem|2]is also true in higher dimensions. This follows 
from the fact that the terms Q/| |z| j'^ for different /3's are orthogonal in the 
sphere measure: 



d-l 



?i • Q2 = VQi, Q2 harmonic, degQi ^ degQ2 
and in particular Q is orthogonal to 1 if dcg Q > 0. 



Theorem 3. The coefficients of the high-order expansion of the harmonic potential a of 
the regular random walk on 1? can he calculated from the differential equations and (|6j. 

Proof. There is almost nothing to prove here. We saw during the proof of the pre- 
vious theorem that | |23> allows to calculate all coefficients except those where /c = 
or / = (i.e. the coefficients of the harmonic summands). There is only one of these 
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(and its conjugate) in every level k + I ~ — n. Therefore the coefficient of z " can 
be determined from the coefficients of z^"~^z, z~"~^z^, . . . and from □ 

The proof of theorem|3]is also the second algorithm for getting l|3). 

Theorem 4. 

(1) The term of order — fc in any walk is P/| |z| |" where a < 4fc + 2 — d and P is a 
polynomial of degree a — k<3k + 2 — d. 

(2) If D3 = 0, in particular if R is reversible^, then the term of order — fc z's P/||z||" 
where a < 3k + 2 — d and P is a polynomial of degree a— k<2k + 2 — d. 

Proof. Take as an example clause 2. We may assume that the walk is spherical. In 
two dimensions, the theorem is equivalent to showing that only terms of the form 
z'^'z', < fc < ^\l\ (or < Z < -11^1) appear. On the one hand, k can only increase 
(by 1) in the step of taking D^^ (see the proof of theorem|2), and this operation also 
increases Z by 1 . On the other hand, the step of taking D4+Dc, + - ■ ■ decreases I by at 
least 4, and the theorem follows. In higher dimensions the same argument works 
using the decomposition into spherical harmonics discussed before theorem|3| □ 

We note that clause 1 of theoremlHis already mentioned in IFU96IIU981 . 

Theorem m can be generalized to cases where higher symmetries exist. For ex- 
ample, for the 6-walk on the triangular lattice, D3 = 0^ = and D4 ~ j'^c which 
allow to get that the term of order |z|^'^' is P/\z\'^ where v < ^k and deg P < 

Theorem 5. The harmonic potential of a reversible random walk contains only terms of 
even order in even dimension and only terms of odd order in odd dimension. 

Proof. Since -D2fe+i = for all k, we only need to handle the terms we cannot 
calculate, namely the (continuously) harmonic terms. However, a harmonic term 
of these orders must be P/||z||" with deg P odd, and since a(z) = a{—z) such terms 
cannot appear in the expansion of a. □ 

4. Explicit constants 

This note is a spin-off from a different work [K] where it was necessary to have 
an explicit value for the constants in the 0( ) in 0. How should one go about to 
calculate something like that? Generally, one would try to get some explicit con- 
stant in the O of an approximation of one order higher, and then use the McCrea- 
Whipple algorithm to check the first few values. 

Which approach would be best to get an explicit constant in the O ( • ) ? The differ- 
ential approach doesn't seem to give any explicit constants whatsoever. Therefore 
it was necessary to use the direct approach. A program was written to handle all 
the algebra and the calculations of the explicit constants in the 0{-)7 After careful 
selection of parameters, the program got ^ with an error term which is smaller 
than < 10^^|z|^^^ for any |z| > 1200. The program doing the algebra ran for half 
an hour on a PC. This may or may not be enough — it depends on the difference 

^ A reversible (or symmetric) walk is a random walk satisfying R= — _R. In graph language, a reversible 
random walk is a walk on a simple graph, while a non-reversible walk is a walk on a directed graph. 
Many of the claims of this paper need slight alterations for a non-reversible walk. For example, we 
have Aa = <5{oj where A is the Laplacian of the time-reversed walk. 

The program is several thousands of lines long. It is available from the first author for examination, 
further development etc. 



14 



GADY KOZMA AND EHUD SCHREIBER 



between the actual maximum and the asymptotic value. However, this difference 
is > 0.01 (see below) and lO^^jzl^^^ <• io-5|2|-2 for \z\ > 1200 so this quantity is 
indeed negligible. 

Interestingly, here the third approach (the one discussed briefly on page |3j 
shows a slight advantage over the other two. The third approach allows to get 
an explicit constant in the O(-) in l|8j with much less effort — a few trivial pro- 
grams to find optimal parameters were all that was necessary. Contrariwise to the 
simplicity of the programs, the estimates are not as good and it was necessary to 
run the McCrea- Whipple algorithm for \z\ < 6000, which takes a few days. 

The maximum of jzp (a(z) — ^ log \z\ — A) was found at z = 3 and so the con- 
stant implicit in the 0( ) in J2j is 

432 + 91og72 + 187 

2 Z 1 _ 153 ~ 0.06882 

TT 

Notice that this is only slightly larger than the asymptotics ^ w 0.05305. Natu- 
rally, explicit constants can be calculated for higher order estimates as well. 

4.1. Other simulations. We add below two additional as. The first is for a "drunk 
king's walk" — a random walk on 1? going at each step to each of the 8 neighbor- 
ing lattice points with equal probability. Here are the first few coefficients: 

4 , , , , ^ / Uz-i 5z8 167z8 

a{z) = — log|z| + A + Rc 



37r ^' ' V97r|2p 907r|z|8 367r|z|i2 2527r|z|" 

35^12 1673z** 15zi2 1925zi6 



547r|2|i8 21607r|z|i6 27r|2|20 2887r|2|24 

so in this case it is no longer true that all the coefficients have a constant sign. An- 
other interesting case is the non-reversible walk on the directed triangular lattice 
going at each step with probability ifoz + l, z + worz + ijJ^ where uj — + ~ 
^/\. Here are the first few coefficients: 

, , %/3, , , , V3^ / 2^ 17z^ 
a{z) = — log|z| + A4 Rc . , . 



TT ' TT \Q\z\'^ 12|z|8 18|2|6 54|zP I35|z 

35zi2 I9z9 7zi5 85z^ 982^2 385zi8 
+ + —rr^ + ,-,0 + ,.0 + 



216|z|i6 54|z|i4 I8|z|20 Il34|z|i2 81|z|i8 324|z|24 

We see that the phenomenon of constant sign is not specific to the regular random 
walk (we checked this up to coefficients of order 0(|z|^28^) note that in this 
case all coefficients are positive while in the 1? case all coefficients are negative. 
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